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ABSTRACT 

In this paper the general problem of a shallow shell with constant 
curvatures is considered. It is assumed that the shell contains an arbi- 
trarily oriented through crack and the material is specially orthotropic. 

The nonsymmetric problem is solved for arbitrary self-equilibrating crack 
surface tractions, which, added to an appropriate solution for an uncracked 
shell, would give the result for a cracked shell under most general loading 
conditions. The problem is reduced to a system of five singular integral 
equations in a set of unknown functions representing relative displacements 
and rotations on the crack surfaces. The stress state around the crack tip 
is asymptotically analyzed and it is shown that the results are identical 
to those obtained from the two-dimensional in-plane and anti-plane elasticity 
solutions. The numerical results are given for a cylindrical shell contain- 
ing an arbitrarily oriented through crack. Some sample results showing the 
effect of the Poisson's ratio and the material orthotropy are also presented. 


1 . Introduction 

Because of their potential applications to the strength and failure 
analysis of such structurally important elements as pressure vessels, pipes, 
and a great variety of aerospace and hydrospace components, in recent past 
the crack problems in shells have attracted considerable attention. Typi- 
cal solutions obtained by using the classical shallow shell theory may be 
found, for example, in [l]-[4]. In a Mode I type of shell problem (that is, 
in a shell for which the geometry and the loading are symmetric with respect 
to the plane of the crack), particularly for membrane loading, the solution 
based on the classical theory seems to be adequate. However, in skewsynmetric 
or nonsymmetric problems, because of the Kirchhoff assumption regarding the 

^*^This work was supported by NSF under the Grant CME 78 09737, NASA-Langley 
under the Grant NGR 39 007 Oil, and by the U.S. Department of Transpor- 
tation under the Contract DOT-RC-82007. 
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transverse shear and the twisting moment, in the classical solution it is 
not possible to separate Mode II and Mode III (i.e., respectively in-plane 
and anti-plane shear) stress states around the crack tips. In this case a 
singularity of the form r"^ in Mode II stress state automatically implies 
r singularity in Mode III. For flat plates such drawbacks of the classi 
cal theory was pointed out in [5] where it was found that the asymptotic 
results obtained from plate bending and two-dimensional elasticity could 
be brought in agreement provided one uses a sixth order plate theory (e.g., 
that of Reissner's [6]). 

In the crack problems for shells even though the membrane and bending 
results are coupled, the asymptotic behavior of the membrane and bending 
stresses around the crack tips should be identical to those given by respec 
tively the plane stress and plate bending solutions. This was shown to be 
the case for the classical shell results (see, for example, the review 
article [7]). Recent studies using a Reissner-type shell theory [8], [9] 
shows that similar agreement is also obtained between shell results and 
those given by the plane elasticity and a sixth order plate bending theory 
[10]-[13]. 

Because of the high likelihood of Mode I type fracture most of the pre 
vious studies of crack problems in shells were on the symmetrically loaded 
structures in which the crack is located in one of the principal planes 
of curvature. The advantage of this crack geometry is that one can always 
formulate the problem for one half of the shell only as a symmetric or an 
antisymmetric problem and reduce the number of unknowns. However, in such 
structural components as pipes and pipe elbows, if, in addition to internal 
pressure and bending the external loads include also torsion, then the 
most likely orientation of the crack initiation and propagation would be 
along a helix rather than a principal plane of curvature. In this case, 
the problem would have no symmetry and all five stress intensity factors 
associated with the five membrane, bending, and transverse shear resultants 
on the crack surfaces would be coupled. Consequently, the related mixed 
boundary value problem would reduce to a system of five pairs of dual inte- 
gral equations or five singular integral equations. 
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In this paper we consider the simplest and yet, from a practical view- 
point, perhaps the most important such problem, namely a cylindrical shell 
containing a through crack along an arbitrary direction with respect to the 
axis of the cylinder. In formulating the problem it is assumed that the 
regular solution of the shell without the crack for the given applied loads 
is obtained and the problem is reduced to a perturbation problem in which 
the self-equilibrating crack surface tractions are the only external loads. 

2. The Basic Equations 

The problem under consideration is described in Fig. 1, As in EH]- 
[13] the material is assumed to be specially orthotropic, that is the elas- 
tic constants defined by 

^11 " ^'^ir''l'^22^ * ^22 " ^ (ct22’''2‘^11^ » 

e 12 ~ ^12^^^12 * '^l^^l ~ '^ 2 /^ 2 . (2.1) 

satisfy the following factorization condition^ ' [12] 


2(l+y^v^v^ 

Defining the following "effective" material constants 

E = *^E-j E2 , V = *'v-|V2 , B = ^ , G = 2{Hy,) * ^ ~ » 


( 2 , 2 ) 


(2.3) 


^ ^The results given in [12] show that the effect of material orthotropy 
on the stress intensity factors can be quite significant. In practice 
the material may be orthotropic because it is either a composite laminate 
or a rolled sheet metal alloy. Orthotropic materials are also anisotropic 
with regard to their resistance to fracture and crack propagation. Hence, 
in a cylindrical shell if the axes of orthotropy do not coincide with 
the axial and circumferential directions, the solution of the general 
inclined crack problem becomes all the more important. The solution is 
also necessary to analyze the weld defects and cracks initiated in the 
weak cleavage plane of the rolled sheet in spirally welded pipes. 
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equations (2.1) may be written as 

1 ^11 1 *^12 
^11 ~ E “ ^^22^ * ^22 ~ E (c2a22“'w^l) . £-12 “ ~2G~ ‘ 

The derivation of the differential equations for a specially orthotropic 
shallow shell based on a transverse shear theory [8], [9] may be found in 
[11]-[13] and will not be repeated in this paper. Referring to Appendix A 
for notation and to [11]-[13] for details, in terms of a stress function (j» 
and the z-component of the displacement w the problem may be formulated as 


follows: 

VW H- X^(1-.V^)(X^ ^ ^ + X^ 

= X-(1-kv2)(^) , (2.6) 

- i{)-w = 0, (2.7) 

72f2 - n = 0 . (2.8) 


The shell parameters X 2 , X 12 , x, and k are defined in Appendix A, 
q(x,y) is the transverse loading, and the curvatures are given by 

J_ = _ _L = . a^Z fo q 

R-J “ 3X^2 » ax^^ * R-j2 " 3 Xi3X2 ’ ' * 

where Z = Z(xi, Xg) is the equation of the middle surface of the shell. 

The functions tp and 0 are related to the components of the rotation vector 
by 


o = Izy. an . _ 3i|) 1-v an 

^x ax 2 ay ’ "^y ay 2 ax * 


(2.10) 
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The normalized membrane, moment, and transverse shear resultants are given 
by 


N = N = N = - • 

XX ’ yy 3)^ ’ xy 3x3y * 


( 2 . 11 ) 


3B 30U 

”xx “ k!^ <1# • V “ iiF 


96, 


3B 

An 


3X * 


M 


_ a 1-v /A + An . 
- - t - > 


36. 


xy ~ hx^ 2 ' 3y 


3X 


V =1^+6 ,V =|^+B . 
X 3x *^x y 3y ^y 


(2.12) 

(2.13) 


3. General Solution of Differential Equations 

Eliminating (j) in (2.5) and (2.6) one obtains an eighth order differen- 
tial equation for w(x,y). If the solution of this differential equation is 
expressed as 

oo 

w(x,y) = f(x,a)e’^^“da , f(x,a) = R(a)e'"^ , (3.1) 

*00 

the characteristic equation for m is found to be 
D(m) = m®-(4a^+KXpm® - 4 kXj2^2“^*’’’'^ 

+ [6a‘*+ic(4Xj2'*'2XjX2+X2)a^+X2]m‘* 

+ 4Xj2[^2^K(Xj+X2)a^]aim3 

- [4a'*+Ka^(4X j 2'*’2X jX 2+X ^ )+4X j2'*’2X ^ X2]o2ni2 

- 4XiXj2(l+Ka^)a^im + a‘»(aHicXia2+Xj) = 0 . (3.2) 

In (3.2), by substituting m=is it may be seen that 
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(3.3) 


D(1s) = z a. (a)s*^ = 0 , 

0 •' 

where the coefficients aj^ are real and, hence, the complex roots are in 

conjugate form. Since Re(mj) = Im(Sj), (j=l,...,8), by ordering the roots 

m. of (3.2) properly it may be shown that they have the following property: 
J 


Re(nij+4) = -Re(mj) , Re(mj) < 0 , (j=l,...,4) . (3.4) 

Considering now the regularity conditions at x=+~ from (3.1) and (3.4) it 
follows that 


f(x,a) 


^ . , '"iX 

E R-(a)e , X > 0 , 

1 

8 , , m^x 

E Rj(a)e , X < 0 . 


(3.5) 


Similarly, if we let 
00 

<|)(x,y) = ^ j g(x,a)e"''^“ da, (3.6) 

<•00 

from (2.5), (3.1) and (3.5) we obtain 


g(x,a) = 


E R,.(a)K,.(a)e , X > 0 

J J 

8 m^x 

E R--(a)K.(o)e , X < 0 

5 J J 


(3.7) 


where 


Kj(a) 



(3.8) 


Also, by assuming that 
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(3,9) 


CD 

n(x,y) = ^ I h(x,a)e"''“^ da , 

•00 

00 

I|»(x,y) = ^ I 0(x,a)e"’“^ da , 

— oo 

from (2.7), (2.8), (3.1) and (3.5) we find 
f A-|(a)e*"^^ , X > 0 

h(x,a) = < 

( A 2 (a)e^^^ , x < 0 , 


= -r2 = , 


e(x,a) = 


(4 Rj(a) mjX 
2 — T 6 

1 icPj-l 

8 ^j(“) mjX 


X > 0 
X < 0 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


The expressions for the stress, moment, and transverse shear resultants 
may be obtained by substituting from the solution given above into (2.10)- 
(2.13). The results are given in Appendix B. 

Since the problem has no symmetry, the preceding analysis would give 
its solution for the half regions x>0 and x<0 separately, and since each 
half would have five boundary conditions at x=0, ten unknown functions are 
needed to account for these conditions. Quite apart from the crack prob- 
lem, Ri,...,Rg, A<j, and Ag are the ten unknowns which may be used to solve 
the two half shell problems. In the crack problem following are the con- 
tinuity and boundary conditions which must be satisfied at x=0: 

^ ’ -»<y«» , (3.14) 

’ -«<y«» , (3.15) 
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Nxy(+0.y) = N^(-0,y) , -co<y<o» , 

(3.16) 

Mjjy(+0,y) = M^y(-0,y) , -»<y<» , 

(3.17) 

Vj^(+0,y) = Vj^(-O.y) , -<»<y«» ; 

(3.18) 

Nxx(+o»y) = F^(y) » \y \<^ » 'I 

u(+0,y) - u(-0,y) = 0 , |y|>*^ , ^ 

(3.19) 

Mxx(+o>y) = |y|<»^ » 'I 

- 3j^(-0,y) = 0 , \y\>r/c , ^ 

(3.20) 

Nxy(+0,y) = FgCy) , |y|<i^ , ^ 

v(+0,y) - v{-0,y) = 0, Iyj>i^ , ^ 

(3.21) 

Mxy(+0,y) = p4(y) > \y\<yc , 'I 

By(+0,y) - 3y(-0,y) = 0 , Iyl>v^ J 

(3.22) 

Vj^(+0,y) = Fg(y) , |yl<v^ , ^ 

w(+0,y) - w(-0,y) = 0, |y|>»^ , J 

(3.23) 


where Fi,...,Fg are the known crack surface loads in the perturbation problem 
under consideration. 

From the expressions of the stress, mogient, and the transverse shear 
resultants given in Appendix B, it may be seen that the homogeneous relations 
(3. 14) -(3. 18) can be used to eliminate five of the ten unknown functions 
A-j, and A 2 . The remaining five may then be obtained from the 
mixed boundary conditions (3.19)-(3.23) . The problem may be reduced to a 
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system of five dual integral equations by obtaining expressions for the 
displacements and rotations similar to that given in Appendix B and then 
by substituting into (3.19)-(3.13) . Hcwever, in the problem under con- 
sideration this procedure would be extremely lengthy. A somewhat more con- 
venient approach is the reduction of the mixed boundary conditions directly 
to a system of integral equations. From the nature of the mixed conditions 
it is clear that the integral equations will be singular. In order to avoid 
strong singularities in the resulting integral equations it is necessary 
that the new unknown functions be selected as the derivatives of the "dis- 
placement" quantities rather than the displacements and rotations. Of the 

"displacements" which appear in the mixed boundary conditions B , B , and 

A y 

w may readily be expressed in terms of Rj and A|^, (j=l,...,8; k=l,2) by 
using the solution given in this section. To find u and v we use the basic 
strain-displacement relations for the shallow shells, namely 


'ij 


" 2 ‘^j.i ^ ^.i 


^3,j ‘^3,i^» 


(3.24) 


From (3.24) it can be shown that 

’“l = 2 !fl2 . ^ . 32Z ^ 

3)^ ax^ " 3Xj " 3Xj 3X^ “ 3)^ 3X^ * 


(3.25) 


By substituting from (2.4), (2.9) and Appendix B into (3,25) we obtain 


3 u _ 

3y " 


+ 


2(l+v)Njy, - J <ly + V J -35T •'y 



(3.26) 


From (3.T9)-(3.23) it is seen that the y-derivatives of the relative crack 
surface displacements and rotations are the natural choice for the new 
unknown functions. However, (3.26) suggests that for the in-plane displace- 
ments u and V such a choice would require very complicated analysis. In 
fact, in the present problem it is not feasible to express the functions 
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R. and Ai. in terms of the new unknowns in the d 

J 4* -• + 

selected as the derivatives of u"^-u and v -v . 
unknown functions as follows: 


- Ir„ I ^ If > 


3e„ 3B„ 

Gp(y) = lim - lim , 
^ x->+0 x-v-0 


G3(y) 


Hm r 3V / 2 ^ 3Wi 
*-3y " ^ X ^ ^ 3y^ 


-11m CIt 


G4(y) 


lim - lim , 
x-v+0 x-»”0 ^ 


G5(y) = lim l^-lim 
^ x-H-0 x->-0 


By using now the solution given in this section 
Appendix B, the auxiliary functions G-|,...,Gg mi 
of R-j,...,Rg, A^, and follows: 
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(3.32) 



(3.33) 


If, ,4 ra^RjCa) 8 ra.R.(a) 


a^Kl l-v 


2 ' '""^ ' [Ai(a)-A2(a)]}e do 


CO 

G3(y) * ^ I ^ ^ KjRj[pj+(l+v)o2] - E KjRj[pj+(l+v)o2]}e do , 


(3.34) 


if” 4 Ri 8 R4 

G4(y) =-2j:| [aM E^) 

j V 


2~^^ o i (r^A^-r2A2)]e"^“^ do , 


(3.35) 


w 

®5<y> = - i 1 ^ - 


Z R.]e" da 

5 J 


(3.36) 


Also, by substituting from Appendix B into (3,14)-(3,18) and inverting 
the Fourier integrals we find 


4 8 

E K.R. - E K.R. = 0 , 
2 J w 2 J J 


(3.37) 


4 P,-+(l“v) o^ 8 P-+(l-v)o^ 

E — -r R. - E -=■*— — T- R- 

J KPj-1 J 5 kPj-1 J 


“ 2 ^ ® “0 , 


(3.38) 


4 8 

E m-K.R. - E m-K.R. = 0 , 

ijjj 5 JJJ 


(3.39) 
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4 

2 n ^1 
1 


8 


m.R. 


ia I - ia i j 

5 KPj-1 


+ f (l-v)(a2+r^)(A,-A2) =0 , 


4 8 ^^-jP-i^-f n 

? irV- - 1 - 1 <<(i-'>)(Ai-''2) = “ • 


1 '"’j 


From (3.32)-(3.41) it then follows that 

4 8 


X2 2 


(A-|-A2)/2 = q2(ci) , 

'* 8 

Z p-K-R* - Z p*K*R« ~ ^o(ct) 9 

4 P-f^T 8 P-1^1 

1 KPj-1 5 KPj-1 

8 


-ia(z R. - E R.) = qc(a) , 

1 J 5 J ^ 


where 




qj(a) = I Gj.(t)e''“^ dt , (j=l,...,5) . 


(3.40) 


(3.41) 


(3.42) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 

(3.47) 


In defining q. by (3.47) it is assumed that G.(y) = 0 for |y| > Vc. From 

J w 

(3.37)-(3.46) the unknown functions R-j,...,Rg, A^, and may be obtained 


as follows: 


Rj(ct) ~ z i • (j“l»»»»»8) y 

k“ 1 

Aj(a) = Z Cj|^(a)q|^(a) , (j=l,2) ; 
k™l 

where Bjj^ and Cjj^ are known functions of a. 


(3.48) 

(3.49) 
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4. The Integral Equations 


The relations to determine the functions 6^,..,,Gg necessary to com- 
plete the solution of the problem are obtained by substituting from (3.47)- 
(3.49) and Appendix B into the mixed boundary conditions (3.19)-(3.23). 

From the definitions of G^,...,Gg as given by (3.27)-(3.31) It Is seen 
that they are related to the derivatives of the crack surface displacements 
and rotations. Thus, In addition to requiring that G.(y) = 0 for |y| > ■/z t 
(j=l,...,5), further conditions must be Imposed on these functions In order 
to Insure the continuity of displacements and rotations In the shell for 
x=0, |y| > /c (see (3.19)-(3.23)). That Is, Gi,...,G5 must be such that 


J W " Uj(-o»y)]dy = 0 » 


(4.1) 


where a>., (j=l,...,5) represents the displacements and rotations u, v, w, 

J 

3„, and b . From (3,27)-(3.31 ) the single-valuedness conditions of the form 

A y 

(4.1) may now be expressed as 

^ X 2 X 2 t 

f [S,(t)+{4^) tG5(t)]dt - (4^) f dt [ G5(y)dy = 0 , (4.2) 

-Jc -Jc -S 



(4.3) 

(4.4) 

(4.5) 

(4.6) 
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With the requirements that G.(y) be zero for |y| > /c and the conditions 

J 

(4.2)~(4.6) be satisfied, the second part of the mixed boundary conditions 
(3.19)-(3.23) relating to the displacements and rotations has thus been 
taken care of. The first part of (3.19)-{3.23) relating to crack surface 
loading would then give the integral equations to determine Gi,...,Gg which 
may be expressed as 


1 f It m..x-iay 

lim - ^ a2 z K.R.e ^ da = F,(y), |y| < , 

x>+0 J 1 J J ' 


lim 

x->+0 


_L a 

2tt hT*" 



p.+(l-v)a2 




m.x 
e ^ 


(4.7) 


- ^ (1-v). 


r,x . „ 

iaA^r^e ' ]e"’“^da = F, 


(y). |y| < ^ 


(4.8) 


lim 

x^^O 


00 



•00 


m. K.R.e 

u J 


m.x-i 

J 


“^da = F3(y). |yl 


< C , 


(4.9) 


lim 

X-J-+0 


1 a(l-v) 
2tt hx** 



KP,--' 


m.x 
e ^ 


+ (l-v)(a2+rJ)A^e ^ ]e'''“^ da = -F^(y), [y] < , 


lim 

x-»-+0 



■‘Pj-’ 


m.x 
e ^ 


- f (l-v)aiA^e ^ ]e"’°^da = Fg(y), jy[ < /c . 


(4.10) 


(4.11) 


By substituting now from (3.47)-(3.49) into (4.7)-(4.11) and by changing 
the order of integrations we obtain a system of integral equations for 
Gp . . . ,Gg of the followi ng form; 
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lim 

X++0 


I z Gj(t)dt I 
-/c -o 


V|^j(x,a)e’^^-y^da = F,^(y). (k=l 


..5), |y|<i^ 

(4.12) 


where Vj^j(x,a), (k,j=l,...,5) are known functions. The dependence of Vj^j 
on a is primarily through r^(a), r 2 (a) and the roots mj(a), (j=l,...,8) of 
the characteristic equation (3.2) and is therefore very complicated. How- 
ever, the functions Vj^. depend on x only through the exponential damping 
terms exp(mjX), (j=l,...,4) and exp(r^x) which simplifies the asymptotic 
analysis of the kernels in (4.12) quite considerably. To examine the singu- 
lar behavior of the kernels given by the inner integrals in (4.12) the 
asymptotic analysis of the functions V|^j(x,a) for large values of |a| is 
needed. First, from (3.2) and (3.12) it can be shown that for large values 


of |a[ the characteristic roots mj and r^ have the following asymptotic 
values: 


Pj P 


mj(a) = -|a|(l + 2 ^ - 0^ + . ..) , (j=l 4) , 

p. p.2 

(a) — I a I (1 + 2(j 2. ” •••)* (j ~ 5,. .., 8 ) , 

r^(a) = -laid + - •••) » 


(4.13) 

(4.14) 

(4.15) 


= i«in + Kd-W ■ ^ 


(4.16) 


where Pj = m^ - a^. Then, observing that the coefficients and 
which appear in the expressions of Rj and Aj (see (3.48) and (3.49)) depend 
on a through m^, (j=l,...,8), and rj^, (k=l,2) only, the asymptotic expansion 
of V|^(x,a) for large values of |a| may be obtained by using (4.13)-(4.14) . 
Consider, for example, the integral equation (4,7). By using the asymptotic 
values found for m.(a) and (3.37), (3.39), and (3.42) it can be shown that 

J 

for large values of |a[ we have 


Z KjRj(a) = - 9i(«) • 


(4.17) 
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By adding and subtracting the asymptotic value to and from the integrand, 

(4.7) may be written as 
00 

lira {- f f C«2 5 + SlanMq da 

j^.v+0 J 1 J J ‘♦I I 

00 ^ 

+ ^f e"l°l^ daf G^(t)e’“^^"y^dt}=F,(y), |y|</c . (4.18) 


By changing the order of integration, evaluating the resulting inner inte- 
gral, and then going to limit, the second integral in (4.18) may be expressed 


as 


lim 

x-H-0 


1 

27 



G-|(t)dt ^ I 
0 


-ax 


sina(t-y)da 


= lim 
x-»-+0 


✓c 


y ^ 

-/c 


(t-y)G^(t) 

(t-y')^+x2 


J_ 

4ir 



G^(t) 

“t=y" 


dt 


(4.19) 


The first integral in (4.18) is uniformly convergent, and hence, the limit 
can be put under the integral sign. By substituting now from (3.47), (3.48) 
and (4.19) into (4.18) we obtain 

G-j(t) 5 

J [ + 2 k^j(y,t)Gj(t)]dt = 4TrF-,(y), \y\<'/c , (4.20) 


where the kernels k^^(y,t) are known functions which are bounded for all 
values y and t in the closed interval [-*^, 

Similarly, the integral equations (4.8)-(4.11) can be reduced to 

1 ,2 ^o(t) 5 h 

J ^ k 2 j(y,t)Gj(t)]dt = 4 tt J F 2 (y), Iy|</c , (4.21) 



k 3 j(y,t)Gj(t)]dt = 4irF3(y), 


|y|<»^ . 


(4.22) 
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(4.23) 




J 2 k^j(y.t)6j(t)]dt = 4tt j F^Cy), Iy|<»^ . 

-1^ 




r Gc(t) 5 

J C I k5j(y,t)Gj(t)]dt = 4ir FgCy), Iy|<»^ . (4.24) 

-i/c 


The expressions of the Fredholm kernels, k.j(y,t), (i,j=l,...,5), are given 
in Appendix C. The details of the analysis may be found in [14]. The system 
of singular integral equations (4.20)-(4.24) must be solved under the addi- 
tional conditions (4.2)-(4.6). They may be solved in a straightforward 
manner by using the Gaussian integration technique (see, for example, [15]). 
The major work in this problem is the evaluation of the Fredholm kernels 
k,-,-(y,t), (i,j=l,...,5) which are given in terms of Fourier integrals. To 
improve the accuracy the asymptotic parts of all integrands are separated 
and the related integrals are evaluated in closed form. The details of this 
analysis may also be found in [14]. 


5. Asymptotic Stress Field Around the Crack Tips, Stress Intensity Factors 

For the numerical solution of the system of singular integral equations 
(4.20)-(4.24) the interval {-/c, is normalized by defining 

T = t/^/c , n = y/^ , 5 = x/.^ , 

H.(t) = G.(TV^,(j = 1,...,5) , -1<T<1 . (5.1) 

J J 

We now observe that the index of the system of singular integral equations 
is +1 and its solution is of the following form; 

Hj.(t) = hj(T)/(l-T2)^ , (-1<T<1) , (j=l 5) , (5.2) 

where h-|,...,hg are unknown bounded functions. The membrane, bending, and 
transverse shear resultants may be obtained by substituting from (5.1), 
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(5.2), and (4.47)-(4.49) into the expressions given in Appendix B. The 
asymptotic behavior of the stress field around the crack tip could then be 
obtained by using the asymptotic expansions of mj and rj^, (j=l,...,8; 
k=l,2) given in (4,13)-(4.16) and the following asymptotic relation [16] 

j ^ e’- d. = (h(1)exp[1(a - ^ )] 


+ h(-l)exp[-i(a - 4 ~ | ~ j ) ] + » 

giving 

ho ( 1 ) f” I _ 1 

= 5 [ e““'^' cos[o(l-n) - ^]da 


h,(l) 
4v^2tT ^ 

00 

sinWI-n) - f]da. 

hod) , 

0 

r” p-“UI 


- — 

6 ’ * 

cos[a(l-n) - ^]da 

Zt/Z^ ^ 

0 

hid) 

09 

[ — (l-a| 

[5le““l^l sin[a(l-Ti) - J]do , 

4/2ir * 

0 


hjd) 1 

00 

' — (l-a| 

|cl)e““I^I sin[a(l-n) - f]da 

41'^tT J 

v^T 

0 

"i(’) 1 

<o 

’ ^ 1 

CosCa(l-n) - |-]da , 

5 v7e-“l 

4/27 J 

0 


(5.3) 


(5.4) 


(5.5) 


(5.6) 
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. h ^"4 


5h/i(’) 


4y?7T 


>■“1^1 cosCa(i-n) - |0<lo 


J!_!^ 

12a . >K— 




00 

f -L (1-K.|5|)e-“|5| slnwi-n) - 50da. 


”yy“W^l 


h ( 1 ) 

f ^ sin[a(l-n) - |-]da. 


4v^2Tr ^ VcT 

0 


»xy = I 7= 0-»lc|)e-l«ls1nWl-n) - |]d„ 


4»'^ ^ yfa 


. 1 h M') 

^ 4.^ 


00 

j ^■|£|e•“l5lcos[c.(^-^) - f]do . 


, h5<’) r e-“l5l . 


V = — 

^ 2.^ 


sin[a(l-n) - |-]da , 


hcO) ^-ctkl 

V = - ■ - — - f cos[a(l-n) - T]da . 


Defining the new coordinates r, 9 in n» € plane by 

5 = r sine , n-1 = r cose , 
and using the relation [16] 

fV’ e- £=j;Krz)dz = („ Un'’ f) 

(s>0 t y^O) } 


(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 
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equations (5.4)-(5.11) can be reduced to the following form 

1 e . 1 50 


- hgO) 

N = — 

4^ 


L- ^ sin ^ ^ sin - 


” 1 '” r 5 e 

L 4 COS 2 


4>^ 


Icos 


yy 


^3^^^ r 7 0 1 r3 0 4. 1 501 

L- ^ sin ^ j sin- 5 -J- [y cos y + 7 cos - 5 -J, 




4/^ 


2 4 


f^xy=- 


^’xx=- 


^3^^^ r3 rnc 0 4. 1 rnc 50T 1 ' 0 4 . 1 c-in 

4 ^ ^4 2 + 4 COS 2 > sin 2 4 2^» 


5e- 


hr 1 ,.4n 0 + 1 50-, ^2^^^ 

W" 4 2"^4 


4^?r 


4»^ 


h r5„__0 1 ^^^ 501 

12a^4^°^2 " 4^°^ 2 ^ 


M = 

yy 


M = 

xy 


''x' 


'’ 4 ^'' h r 7.^.9 l.^_59n M’* h r3 0 , 1 56, 

■ w' 5^ ^ w4">=2 

h r3_,,9 . 1.,,56i '’z*'' hr Ir-i-e I Ur„5ei 

■ 125^?'“=? * Tzit- 5S’"z + 4Si»r^> 

•’5<1) 9 

— = cos %■ 

2/Ir 


2 * 


.. - . 0 
V„ = sin -K- . 

y 2v^ 2 


By observing that the membrane and bending components of the stresses 
given by (see Appendix A) 


m _ M u _ icai M /i 4 _„ „\ 

"^ij “ ^■o, ij ■ h ”ij 


b _ 12az 


from equations (5.14)-(5.19) for the leading terms of the combined in- 

h 

plane stresses a- - = a--"’ + cr.. , (i,j=x,y) one obtains 

■ w • J 1 0 


gJ » 

(5.14) 

(5.15) 

(5.16) 

3 

(5.17) 

(5.18) 

(5.19) 

(5.20) 

(5.21) 

are 

(5.22) 


- 20 - 



h3(l)+zh4(l) 


XX 


4v^ 


[- Is1n| + ls1n^] 


h^0)+zhg(1) j. 5 0 1 cos ^ 1 

4/2F C 4 cos 2 - 4 cos 2 ] » 


yy 


■ 4^ f- 4 2 • 4 2 ^ 

h^(1)+zh2(1) p 3 £ 1 50 -, 


xy 


h3(l)+zh4(l) p 3 0 1 50 -, 

— [ ^ cos ^ + 4 cos ^ j 


4»^ 

hT(1)+zh2(1) 

4,^ 


r 1 ® J. ^ e,-{n 50 1 

[- 4 sin 2 4 2 ^ 


Similarly, for the transverse shear stresses from 


'iz “ I ''iC’ - ’ 


we obtain 


••xz 


3 ^5^^^ 0 n / az 

2 2^ 2 ■ ^h/2^ ^ * 


hs(l) 


3 "5'-' 0 


az ^2. 

'yz 2 2^ 2 ■-■ - W 


sin 7 [1 - (tttt) ] • 


(5.23) 

(5.24) 

(5.25) 

(5.26) 

(5.27) 

(5.28) 


Note that for the isotropic materials c=l and the asymptotic stress 
fields (5.23)-(5.25) and (5.27)-(5.28) found from the shell solution are 
identical to those given by respectively the in-plane and the anti-plane 
elasticity solution of a two-dimensional crack problem. If we now define 
the Modes I, II, and III stress intensity factors (for a crack along x-j=0, 
-a<X 2 <a) by 
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(5.29) 


“ lliD i^ 2(^X2~3 ) cTi j (O »X2 » (j"! »2»3) ^ 


X2->a 


from (5.23)-(5.25) and (5.27) and Appendix A we obtain 


~ ^ ^2(1)] » 

k2(x3) = - ^Ch3(l) +^h4(l)] . 

“ “ T bp(l)[l - (h^) 3 . 




(5.30) 

(5.31) 

(5.32) 


6. The Results and Discussion 

The main interest in this study is in evaluating the stress intensity 
factors in shells for various crack geanetries and loading conditions. For 
each crack geometry the problem is solved by assuming only one of the five 
possible crack surface loadings to be nonzero at a time. For a general 
loading the result may then be obtained by superposition. From (5,30) and 
(5.31) it is seen that the in-plane stress intensity factors k^ and k2 
have a "membrane" and a "bending" component, and h-j and hj are related to 
the membrane and h2 and h4 are related to the bending stresses. For sim- 
plicity, the related stress intensity factors are defined separately. The 
calculated results are normalized with respect to a standard stress inten- 
sity factor a.v^ where a. stands for any one of the following five nominal 

J \J 

("membrane", "bending", in-plane "shear", "twisting", and "transverse shear") 
stresses: 

= N^^/h , ajj = 6M^^/h2 , , 

= 6M^2/h2 , = (3/2)V^/h , (6.1 a-e) 

where crack lies in X2X3 plane and Mp, N-J2* ^-^ 2 * ^1 measure 

or amplitude of) the crack surface tractions. 
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The normalized stress intensity factors are then defined and calculated 
in terms of h^.(l), (i=l,.,.,5) as follows: 


■Siu 


_ n 


k,(0) 




"J 


( 6 . 2 ) 


_ k^(h/2)-k^(0) _ cE h . 

• 


kp(o) r 




_ k2(h/2)-k2(0) _ E h . /,x 
k*4 - - - 2a * 

a.^ j 


'tj 


kgCO) 

a . 

w 


= “f~v^h 5(l) , (j=m,b,s,t,v) 
J 


(6.3) 


(6.4) 


(6.5) 


( 6 . 6 ) 


where for each individual loading o. is given by (6.1). In the case of 

w 

uniform crack surface loads N^-|, N^2» ^12* '^l* ’^e^'erring to (3.19) 

(3.23), Appendix A, and (6.1) the input functions of the system of integral 
equations (4.20)-(4.24) are given by 


" cE ’ ^2^^^ " ^ * ^4^^^ " 6E * 

F5(y)=|-^. (6.7) 


Even though the formulation given in this paper is valid for any shell 
with constant curvatures l/R-i* 1/R2» and 1/R-J2* the results are obtained for 
the practical problem of a cylindrical shell containing an arbitrarily 
oriented crack only (Fig. 1). The crack is assumed to be in a plane defined 
by the angle 3 shown in Fig. 1. For the shallow cylindrical shell the 
curvatures referred to x-j, X2 axes shown in the figure and defined by (2.9) 
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Table 1. Stress intensity factor ratios in an isotropic 

cylindrical shell containing an inclined crack under 
uniform membrane loading N-|i, v = 0.3, B = 45®. 
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may be expressed as 


1 _ sln^e 1 _ cos2$ 1 _ sin0cos0 m m 

^ " R * R2 R“ * R^2 ' ' R 

Some numerical results obtained for an isotropic cylinder are shown 
in Figures 2-11. Figures 2-6 show the primary stress intensity factor ratios 

*^nm* *^bb* *^ss* *^tt* *^vv ® cylinder having a crack inclinced 45® 
with respect to the axis. The unusual results here are those found for 
and k^^. Under a twisting moment H-jg uniformly distributed along the 
crack, the Mode II stress intensity factor ratio k^^ appears to be nearly 
independent of the shell curvature 1/R but highly dependent on a/h. Fig. 6 
shows that the monotonic variation of the stress intensity factor ratios 
with a/h and h/R observed in Figures 2-5 and in previous shell solutions 
is not valid for k^^. This seems to be the case for all values of 0 
varying from zero to ninety degrees. 

The effect of 0 on the primary stress intensity ratios k^^, k^^, 

ktt» and k^^ is shown in Figures 7-11. Extensive results giving all stress 
intensity ratios k.. (i,j=m,b,s,t,v) for 0=0, 15°, 30°, 45°, 60°, 75°, 

• J 

90° and for varying h/R and a/h may be found in [14]. Table 1 shows some 
sample results regarding the secondary stress intensity ratios in a 
cylinder with a 45° crack under torsion (i.e., N-j-j = constant and all other 
crack surface tractions zero). 

The stress intensity factors given in Figures 2-11 and in Table 1 are 
obtained for the Poisson's ratio v = 0.3. Some sample results showing 
the effect of v on the stress intensity factors are given in Table 2. 

It is seen that this effect is not really significant. 

It should be noted that the Poisson's ratio v in isotropic shells and 
V = /viV2 and the stiffness ratio c = (E^/E2)^ in specially orthotropic 
shells appear in the expressions of the kernels of the integral equations. 
Thus, to investigate the effect of the material orthotropy on the stress 
intensity factors both v and c must be varied. However, as seen from 
Table 2 the influence of v is rather insignificant. Therefore, to study 
the effect of the material orthotropy it may be sufficient to vary c only. 
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Table 2. The effect of Poisson's ratio on the stress inten- 
sity factor ratios in an isotropic cylindrical shell 
containing an inclined crack, B=45®, a/h=2, h/R=l/10. 



0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

*^rnni 

1.166 

1.167 

1.167 

1.167 

1.166 

1.164 

*^bni 

0.063 

0.077 

0.092 

0.108 

0.124 


'‘sm 



-0.059 




*^tm 

-0.009 

msm 

HEH 





BBa 



-0.070 



*Snb 

0.018 

0.023 

0.028 

0.034 

0.039 


•^bb 

mmm 

0.617 

0.626 

0.632 

0.634 

0.631 

■^sb 

ssm 

msm 


mygn 


-0.026 

•^tb 

mggn 

mmm 



-0.006 


^b 

0.004 

0.003 

0.003 

0.003 

0.003 


*Sns 

-0.058 




mi^n 


•^bs 


BBiia 





*^ss 


1.059 

1.059 

1.059 

1.058 


•^ts 

0.007 

0.007 

0.008 

0.008 

0.009 


*^vs 

0.133 

0.131 

0.129 

0.128 

0.126 

0.124 

•Srit 

0.005 

0.005 

0.005 

0.004 

0.004 


wm 


mmm 









-0.006 

BiBiiiOi 


•^tt 

0.309 

0.325 

0.339 

0.353 

0.366 

0.379 

•^vt 

-0.095 

migiii] 




-0.088 

'Sl.V 

-0.223 


-0.266 




mm 


0.001 

0.007 

0.014 

0.022 

IjBjgJ 

*^sv 

-0.174 


mgin 

Msm 


-0.238 

•^tv 

1.138 

1.166 

1.191 

1.213 

1.233 


Kvv 

2.304 

2.287 

2.272 

2.258 

2.244 

2.231 


- 26 - 






































































































Table 3. The effect of material orthotropy on the stress 
intensity factor ratios in a cylindrical shell 
containing an inclined crack; 6=45, a/h=3, 
h/r=l/10. 
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For a strongly orthotropic material (graphite-epoxy composite) this effect 
is shown in Table 3. The axes of material orthotropy are along 45® direc- 
tions with respect to the cylinder axis and the crack is located along one 
or the other axis of orthotropy. The Poisson's ratio is v = /v-jV 2 = 0.037 
for the orthotropic shells and v = 0.3 for the isotropic results included 
for the purpose of comparison. The table shows that the effect of material 
orthotropy on the stress intensity factors could be very significant. 

The quantity which is of some interest in certain fracture studies is 
the rate of internally released or externally added energy per unit fracture 
area created as a result of-crack propagation. If U is the work of the 
external loads, V the total strain energy, and A the fracture surface, 
then in a quasi static problem the rate of energy available for fracture 
would be d(U-V)/dA. For elastic problems this energy rate is known to be 
the same for "fixed grip" and "fixed load" conditions. It can therefore be 
calculated as the crack closure energy under fixed grip conditions. Under 
these conditions, dU = 0 and for a crack going from Xg=a to X 2 =a+da, dV may 
be expressed as 

h/2 da 

dV = - 


f 1 3 

J I cTij(0,X2,X2)[uj(+0,X2“da,X3) - Uj(-0,X2-da,X3)]dx2dx3 


where the minus sign is due to the fact that during the "release" of the 
crack surfaces in a<X 2 <a+da, -h/ 2 <X 3 <h /2 the tractions and displacements 
are in opposite directions (consequently, the total strain energy of the 
shell decrease). For small values of da we now observe that 

k.(x3) 

a,.(0,Xp,X3) = ■ , (j=l,2,3) (6.10) 

^ ^ /2Ix2-a) 

Uj(+0,X2“da,X3)-Uj(-0,X2“da,X3) = — >^2ta+da-X2) » (j=l,2) (6.11) 

_____ 

U3(+0,X2-da,X3)-U3(-0,X2“da,X3) = — /2(a+da-X2) » (6.12) 


- 28 - 



where k-j, k 2 » and k 2 are the Modes I, II, and III stress intensity factors 
around the crack border X 2 =a, 

Referring to the definitions of the stress intensity ratios k^.j 
(i»j=m,b,s,t,v) , given by (6.2)-(6.6) we can define the "membrane", 
"bending", "shear", "twisting" and the "transverse shear" components of 
the stress intensity factors at the crack tip X 2 =a as follows: 


k, = z k..a.>^ , (i,j = m,b,s,t,v) , 
1 j iJ J 


(6.13) 


From (5.30)-{5.32) , (6.2)-(6,6), and (6.13) the stress intensity factors 
may then be expressed as 


'^1^^3^ = '^m '^b W ’ 

X 3 

k2(^3) *^s *^t ^h7I^ ’ 

k3(x3) = kyCl - • 


(6.14) 

(6.15) 

(6.16) 


By substituting from (6.10)-(6.12) and (6.14)-(6.16) into (6.9) we obtain 

. 2 


k 2 


dV 


ir n. 2 , b , I. 2 , t , 4(l+v) I, 

= “ E L'^m ^ ''s ^ 


(6.17) 


Observing that hda = dA, for the rate of externally added or internally 
released energy (at one crack tip X 2 =a, per unit shell thickness, per unit 
crack extension in the plane of the original crack) we find 


d 

iA 


(U-V) = f [k^2 + 



k 2 


k 2 


+ 


4(l+v) 

15 



(6.18) 


Finally it is again worthwhile to remember that alj[^ shell theories 
are, to varying degrees, approximations of the three dimensional elasticity. 
Therefore, even if the "shallowness" assumption is satisfied, the theory 
used in this paper and the results given are only approximate. Strictly 
speaking, the crack problems considered in plates and shells are three- 
dimensional elasticity problems. Such problems in their simplest form do 
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not seem to be as yet analytically tractable. However, from a structural 
viewpoint, the shell solutions can be useful in the sense that the "plane 
stress" crack solutions are, that is, the results should be interpreted and 
used in a certain thickness-average sense. Since the shell theories are 
quite numerous, there is always the question as to what theory to use in 
the crack problem. Clearly there is no unique answer for this question. 
However, one could try to establish some guidelines and set certain mini- 
mum requirements. In crack problems the most important information (from 
an application viewpoint) is imbedded in the asymptotic solution of the 
problem around the crack tips. The first requirement then is that the 
asymptotic results found from the shell solution must be compatible with 
that of the in-plane and anti-plane elasticity solutions of the crack problem. 
This means that the stresses around the crack tips must have the standard 
square root singularity and their angular distribution must be identical 
to that given by the related two-dimensional elasticity solutions. 

In crack problems since one is interested in the behavior of the solu- 
tion very near the crack tip, it is natural to assume that all local length 
parameters would have some influence on the results which are of interest. 

In a general shallow shell there are five local length parameters, namely 
three radii of curvature, R^, Rg, R-| 2 » the crack length 2a, and the thick- 
ness h. A particular shell theory to be suitable for crack problems should 
therefore contain four dimensionless (independent) length parameters. 

Again, since it is desired that the shell theory give a reasonably 
accurate solution near the crack tip, it would be necessary that the theory 
should accommodate all the stress boundary conditions on the crack surfaces 
separately. 

Reissner's transverse shear theory, which has been used in this paper, 
seems to be the simplest theory which satisfy all these requirements. 

Aside from a certain degree of confidence one may have in its results, an 
advantage of such a compatible theory, is that it makes it possible to carry 
out calculations such as that of energy release rate (see (6.18)) routinely. 
This, of course, is primarily due to the fact that the asymptotic results 
(5.23)-(5.25) and (5.27)-(5.28) are identical to that of the corresponding 
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elasticity solutions. However, since a higher order shell theory does 
not necessarily imply higher accuracy in (certain calculated) results, 
there are still unresolved questions. Are the results of the crack problems 
obtained from the Reissner's shell theory, for example, more reliable than 
that given by the classical shell theory? For the stress intensity factors 
we think the answer is yes. The reason for this is largely the fact that 
the classical theory satisfies none of the requirements listed above. Could 
one improve the solution further by considering "higher order" theories 
which may take into account additional features of deformations and stresses 
(such as, for example, the stretch in thickness direction)? Even if one 
can solve such problems with the same degree of numerical accuracy as the 
problems based on simpler shell theories, it would be difficult to know 
which solution is more reliable. In our view, therefore, it would be very 
difficult to justify the use of a more complex theory than Reissner's in 
solving the crack problem in shells unless one has a denonstrable reason 
for it. 
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Appendix A 


Dimensionless and normalized quantities used in the analysis, 

X = x^/a*^ , y = Xgt^/a , z = X3/a ; 
u = u^*^/a , V = Ug/av^ , w = u^/a ; 

Bx = » 6y = ’ 

<l>(x,y) = F(x^,X 2 )/Eha 2 ; 

_ ^ ^ _ ‘='13 _ '^ ‘='23 . 

*^xx " Ec * ‘^yy T“ * ''xy " E * °xz “ ’ °^yz " B ’ 

Nxx = N^^/Ehc , Nyy = cNgg/Eh , N^y = *• 

\x = M„/Ech2 . Myy = cMgg/Eh2 , Mxy = '* 

= V:,/Bhy^ , Vy = Vgv^/Bh ; 


Xj = 120-V2) = 12(l-v2) ^ 


= 12(l-v2) , x** = 12(l-v2) = E/BX** ; 

E, 

E ~ /E^ ^^2 ’ ~ i>2 > ~ ^ 2 /^^ * ~ E * ^ 


5E 

12(l+v] 


(A,l) 
(A. 2) 
(A,3) 
(A.4) 

(A,5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 

.(A.IO) 
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Appendix B 


The stress, moment, and transverse shear resultants. 




Mxx(x.y) 


1 ^ ni*x i 

- ^ J a2 E KjRj(a)e ^ e"^“^da X > 0, 


—00 

-00 


T f~ 8 4 w 

- I «2 E K.R^(a)e ^ e"^“^da x < 0, 
dTJ } ._C 3 3 


(B.l) 


j=5 

r ^ (i-v)p^+Pj 


, 5 R („)e je- 1 «yd„ 

Ztt hA^ J . 1 '^Pi-1 J 

-00 ^ ^ 

- |ifl,(»)c.r,e''’V’“)'da x> 0 . 


8 (l-v)a2+p. m^x 

r aL 

=5 






x<0 . (B.2) 


Nj^(x,y) = 


“ji/j 


m.x . ,, 

•m.e e"^“'^dct x>0 , 

w w 


i f a E R^(o)K^m.e ^ e"’*“^da x 0. (B.3) 

J j=5 J J J 
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Mj^y(x,y) = 


Vj^(x,y) = 




” 'h ^ e“"‘“^do. 

*00 

x> 0 , 

1 n-v)a fT 2 _"'j*.-iayj 

-ZJ-hix^ .fcipTT-® ® 

*00 ^ ^ 


x<0 . (B.4) 


1 r 4 KP,-R,-(a)m- m.x _. 

- e''’V^“yda, x>0 . 


j.r ? 

L j=5 ■'Pj- 


S ^Ve-i»y d. 


[ - j- I iAjCajo e'^^V’“^da^<0 (B.5) 
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Nyy ( X . y ) = 


Myy(X,y) = 


Vy(x,y) = 


1 f 4 

^ r R. tCra.^e'^ da, x> 0 , 

27 t J . ^ 3 3 3 

-00 ” 

^ j .fg x<0 . 

_M — 


(B.6) 


~ 4 vP.-(l-v)a 2 «X 


J_ a r J 

I 2 tt hA^ ^Pj"' 

57 A,(a)ar,e' 


. e e~^“^do 


r^^x g-iayj^ , x> 0 . 


r 8 vPj-(1-v)g^ nijX ^„y 

A ■■’Tp^ “j® ® 

_**» J ^ J 


'*’ jiT F iA 2 (ct)ctr 2 e ^ e"’“^da, x<0. (B.7) 


da 


f - -X — I ia Z — ^ r R- G 
~^l A^(»)r^e ' e”’“^dg, x>0 , 

— oo 

1 8 •«?,• .iav 

■ 57 J A 

_« J 

[ - ^ f A 2 (g)r 2 e'' 2 V’“y da , x<0 . (B .8) 
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APPENDIX C 


The Kernels of Integral Equations; 


CO 

T,»^ n) = I [2a^Re( E KjBj^-l]sincti^(T-Ti)clo 

f" 8 

+ 2 I a^lfu KjBj-j)cosct»^(T“n)dct , 

r ® 

^12^*^ T,»^ n) = 2 I a^ReC E KjBj^sinai/c(T-n)doi 
f“ ® 

+ 2 J a^Im[ E KjBjJcosay^(T-n)da » 


00 

kisCi^ Tf/c n) = 2 I a^Re[ E KjBj^sina)^(T-n)da 
00 

+ 2 I a^Im[ E Kj^ 2 [lcoS“’^('r"n)da » 

0 ^ ^ 

00 

ki4(»^ T,t^ n) = 2 I a^Re[ E Kj^4]sina>^(T-Ti)do 

oo 

+ 2 I a^Im[ E Kj^ 4 ]cosav^(T-n)da » 


(C.l) 


(C. 2 ) 


(C. 3 ) 


(C. 4 ) 
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n) “ 



o 


T,v^ n) = 



0 


^ 22 ^'^ T,*^ n) = 



0 


k 22 (»^ T,i^ n) = 



o 


{ 


f “ Q 

2 I a^Re[ I K^g]sino»^(T-n)clo 


i^ImC E 
0=1 


Kj ^ 5 ] cosa *^( T - n)da 


(C.5) 



Pi +( l - v ) a ^ 

— Bj ^] sfnoi ^ CT - n)cla 


8 p ,-+( l - v ) a 2 

Ini[ z — 1 B.T3cosai'*c(T-n)do » (C.6) 

j=l ’" Pj "' 


1 


2 r ® P-j+(l-v)a2_ 2kO-v)2ar-j , 2 

^ - ^-1 - V -> 


X sino >^ CT - n)da 
^ 8 P ^-+( l “ v ) a 2 

[m[ z — j^pTiT Bj2]cosa»^(T“Ti)da , (C.7). 

0 "l 0 


r 8 P^+(l-v)a2 

J Bj2]sinai/c(T-n)da 


8 Pj + 0 - v ) a ^ 

fmC z 


I P V] — Bj2]cOSa»^(T-T|)da , (C.8) 

0"1 j 
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^24^*^ n) 

2 

:r 

^25^^^ n) 



kgiCv^ T,i/c n) 


-I 


l<22(*^ T,»^ n) 


+ 2 


r“ 8 Pi+0-v)a2 

J ^ ~^ Tp ~ ^“ Bj JsinaV^CTTi)clc 

_ j“1 : j 


f 


8 P»+(l^v)a2 

, z 1 

j«l "Pj* 


Im[ Z ■ ^ J, ' , 1 Bj^]cosav'c(Trn)da « (C.9) 


2 f” 8 Pi+(l"v)«^ 

” " J ~ ^"kP -1 ^j^sfna»^(T-n)do 


’ 8 Pj+(1-v)o2 

Im[ z - ^ — — ^^cosa»^(T-n)da» (C.IO) 

j“l j 


00 

= -2 I a KjnijBjj]cosa/E'(T-n)d 



8 

a Im[ z K.m,3^-nsinai^(T-n)da 
j=l ^ 


(c.n) 


00 • ’ 

= -2 I a Re[z 


® g 

a In)[z Kjnij^ 2 ]s 1 no/c(T-n)da 


(C.12) 
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co> 

f 8 

1^33(1^ T,j^ n) = -2 I o Re[ L Rjinj^jlcosav'cCx-nJcIa 

^ jal 

+ 1 {2alm[ E Kjinj^3]-l}s1na>^(T-n)cla » (C .13) 

^34^*^ T,/c n) = -2 |"a Re[ E Kjin^43cosaV^(T-n)da 
00 

f 8 

+ 2 I a Im[ E K.mAiilsincti^(T-n)da , (C. 14) 

j 4--1 j jj** 

0 

00 

f 8 

^35^*^ T,»^ n) = ”2 J a Re[ E Kjtnj^glcosa»^(T-Ti)da 

00 

+ 2 I o Im[ E Kjnij^]sina/c(T-n)do » ( C.15) 

k^.j()/c T,i/c n) = j a Re[ E icp^^-l ^3cosav^(T-n)da 


~ — Q - ^ -^ o Im[ E — r^— rBi-i ]sinav^(T-n)da 
X J p.-l J1 


j=l 


(C.16) 
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n) 




a m 


a Re[ E ' TB.,]cosa»^CT~n)da 
j=l 


2(l-v) 


f 8 

J “ ^ ^^p*^Bj2]sinov^(T-n)da , 

o j 


(C.17) 


CD ^ 

^43^*^ T,*^ n) = I ® RsC E ^g3cosai^CT-n)cla 

• n J 


20-v) 

A** 


CO 

0 


O -•••J 

Im[ E — e*— TrB;o]s1nai^(T-n)do , 

j=l •'Pj’* 


C C.18) 


T,*^ n) = j a Re[ E B j^3cosav^CT-n)da 


f- ,2(l-v) T_r 2 ***j R Ti 1“''^ + l+a^icd-v) / 9i„ 2^ 

- J ‘"V"^ “ -P“ * — a>i> 4-- ' - (c^+ri ) 


x(l-v) }sinai/^(T-n)da > ( C.19) 


*^45^*^ T,i/c n) = 1 a Re[ tc ' p’^!-1 ^5]cosa»/c(T-ri)da 

CO 

- ^ I a ;^p^^5]stna/c(x-n)da 

n V ^ f” 

+ \4 * — = cosoi^(T-n)do C C,20) 

^ Jo 
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n) = -j Re[ s ]sina»^CT-n)cla 

J — 1 J 


■f 


Im[ s - 7 fM-EL,]cosai^CT-n)da , 
j=i 


*^52^*^ T,>^ n) = -j Re[.2^ 


j=l ■'•'z 


f~ 8 KP,-ni^ 

-] I m[ s 
- j“l j 


f“> 8 •^Pi^'^-i 

kgaCv^ n) = -J Re[ r ^^^^^ 3 ]sina/c(T-n)da 


j=l 


r 8 x^Pi"ii _ 

-J Im[ s ■j^p'|^^ 3 ]cosa>^(T-n)da 

n J“^ j 


f“ 8 KPi'ni,' _ 

kg4(y^ T,^ n) = -J Re[ s ^^p=y-^4]sinai/c(T-n)da 

_ J ~ ^ J 


[- 8 ^ 

-J Im[ S ^^^|y-^4]C0Sai/c(T-n)da 

o j 

-f” ItSS^Elbi coSa.^CT-nlda , 


(C.21) 


CC .22) 


( C .23) 


CC ,24) 
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kgsC*^ x,»^ n) 


= -f {Re[ z 

i j=l 


KPjm. 

^ 1 )}stnaVcCT-n)dc( 


tcp 



Kp.m. 

^^p7::^j5]cosav^(T-n)da , 


( C.25) 


In the expressions given above Bjj^(a), (j=l,...,8; k=l,,,.,5) 
are the coefficients given in (3,48) which are obtained by solving 
the linear algebraic equations (3,37)-(3.46) for the unknowns 
y • • • >Rq* A"] » and A 2 » 
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R 

Figure 2. Stress intensity factor ratio in an 
isotropic cylindrical shell containing an 
inclined crack under uniform membrane 
loading Nn; 8=45°, v=0.3. 
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Figure 3. Stress intensity factor ratio kbj, in an 
isotropic cylindrical shell containing an 
inclined crack under uniform bending 
moment M-ji; s=45®, v=0.3. 




Figure 4. Stress intensity factor ratio kgs in an 
isotropic cylindrical shell containing an 
inclined crack under uniform in-plane 
shear loading N-jpi 3=45°, v=Q.3. 






Figure 5, Stress intensity factor ratio in an 
isotropic cylindrical shell containing an 
inclined crack under uniform twisting 
moment M-ip; 6=45®, v=0.3. 



Figure 6. Stress intensity factor ratio in an 
isotropic cylindrical shell containing an 
inclined crack under uniform transverse 
shear loading V-j; b=45", v=0.3. 
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Figure 7. Stress intensity factor ratio knjn.in an 
isotropic cylindrical shell containing an 
inclined crack under uniform membrane 
loading Ni]; v=0.3, h/R=l/5. 
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Figure 8. Stress Intensity factor ratio kk^ in an 
isotropic cylindrical shell containing an 
inclined crack under uniform bending 
moment Mn; v=0.3, h/R=l/5. 
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Figure 10. Stress intensity factor ratio k+t in an 
isotropic cylindrical shell containing an 
inclined crack under uniform twisting 
moment M-|2; v=0.3, h/R=l/5. 
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Figure 11. Stress intensity factor ratio kyy in an 
isotropic cylindrical shell containing an 
inclined crack under uniform transverse 
shear loading V-j; v=0.3, h/R=l/5. 
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